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Abstract 

that non topological solitons, Q balls, evaporate into fermions. 
All the constructions we used to find evaporation rate were dased on the fact that no fermions 



f— ^ \ We know from previous work 



^ , would move towards the Q ball. All these constructions left an opened question that is : what 

0^ \ happens when a fermion interacts with a Q ball. We shall answer this question in this work by 

using the constructions done to compute evaporation rates. We shall also obtain a new approach 

to compute evaporation rates. 



I. INTRODUCTION 

A scalar field theory with an unbroken continuous global symmetry admits a remarkable 
class of solutions, non-topological solitons or Q-Balls. These solutions are spherically sym- 
metric non-dissipative solutions to the classical field equations [jj, |2, 14 lll| • These solutions 
are constructed using the fact that they are absolute minima of the energy, this property 
enshures the stability of Q-balls versus decay into scalars. The addition of a coupling to 
fermions will result in Q-balls becoming unstable versus decay into fermions J2; 14 III I- The 
result of this instability is Q-ball evaporation. 

Q-ball evaporation admits an absolute upper bound resulting from the fact that all the 
energy ranges are finate. The fact that energy ranges are finate comes from the mixing of 
positive and negative frequency terms. We could ask, what would happen if the fermions 
interacting with the Q-ball have an energy lying outside the evaporation range? This is the 
question we propose to answer in this work. We know from previous work i4j] that we can use 
two methods to obtain a solution to the equations of motion describing a fermion interacting 
with a Q-ball. Our starting point will be the solutions expressed in terms of partial waves, 
that we shall study outside the evaporation range. This partial wave superposition has 
the reflection and transmission amplitudes as expansion coefficients. So we shall study 
the diffusion of a fermion on a Q ball. We shall then apply the same calculations to the 
total energy spectrum to see to the interference of both phemomenoms (evaporation and 
diffusion) . 

In fact we expect the difficulties to appear when the energy of the interacting fermion lies 
inside the evaporation range for this reason we shall start with fermions having their energy 
bigger than the evaporation upper bound. In this situation expect that the only wave having 
a non-trivial behaviour would be one going through the Q ball. The second regime where 
the energy of the interacting fermion lies inside the evaporation range will lead to some new 
results. 

The work will be organised as follows we shall first give a brief overview of the method to 
construct the solutions. We shell then study these solutions outside the evaporation range 
and then we shall try to look at the total energy range. 



II. PRELIMINARIES 

The details of all the following calculations can be found in previous work (2,14; III | and 
we shall only give here a brief overview of the results and methods used. The major idea to 
solve this problem is considering a Q ball localized in space, so the fermionic field equations 
are free ones outside the zone where we have Yukawa coupling (interaction of Q ball and 
fermions). For the two cases we shall consider, massless and massive fermions, we shall need 
to match the solutions outside to the solutions inside the Q ball. The matching coefficients 
will give the reflection and transmission amplitudes we shall use in the problem. 

A. Solutions to the equations of motion, massless case 

Writing down the Lagrangian of a massless fermion having a Yukawa interaction with a 
scalar field gives in one spatial dimension, 

^ferm. = i'ipa^df.'ijj + {g(f)ip'^ip + h.c), (1) 

Equations of motion for the two components of the \l/ field are : 

{ido + id,)iJi - gcpipl = 0, 
{ido - id^)^2 - 9(t>*^i = 0. 



and = 0oe *'^°* in the zone from — / to +/ and zero everywhere else. Using the ansatz : 
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the equations of motion are reduced to the following 2x2 linear system 

0. 




The determinant of the system gives k = ±Ve^ — M^ = ±k^. Solving for the two cases 
k = +/Ce and k = —k^, we obtain the solution inside the Q-Ball: 



^Q = r^M = A e-'"^' + b\ ^' e''^^ (4) 





where M = g(f)Q, g is the couphng constant and 0o the value of the scalar field. The second 
part is the solution when 0o = (outside the Q-Ball) it is, 
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where superscripts L, R indicate the left and right side of the Q-Ball. In order to solve 
Dirac's equation everywhere the solution needs to be continuous in space. Space continuity 

gives at z = —I : 

C^ = Aa.e^^'^^-^)' + 5e-*(^^+^)', 



and at z = +1, 
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These two relations will give the matrix linking the coefficients on the left to those on the 
right. 

B. Solution to the motion equations, massive case. 

Using the same Lagrangian as for the massless case and adding a Dirac mass coupling 
for massive fermions, gives the fermionic Lagrangian : 

C = ^pi^df^^p + g{^^i^<p + h.c.) + Md{x)^^. (6) 

The equations of motion using two component ■?/'-field, in 1 © 1-dimensions are 

{idt + id,)iJi - Me-^^Vs + MDi^2 = 0, 

{idt - id,)ij2 + Me-^^Vi + Moi^i = 0. 
These equations can be solved if we use fields with four degrees of freedom. Solving this 
system will give us the solution inside the Q-Ball, which is the first step. To solve these 
equations of motion we use the following ansatz 

^i = /i(z)e^(-^)* + /2(^)e-(^+'?)*, 
^2 = ^i(z)e^(^-^)* + ^2(-2)e-^('+^)*. 



Due to the separation of time components this ansatz will lead to four equations. Taking, 



fi{z) = Ae'P', f^{z) = Be'^',gi{z) = Ce'P'g*{z) = De'^'. 



and after some re-arrangement of the equations we obtain : 



(9) 



(10) 



-e-/i -92 + MdQi = pfi, 

^-9i - f2 - MdIi =P9u 
-e+/* + 91- Md92 = Pf2, 

e+(7* + /i + Mof^ = P92, 

where e_ = e — ^ and e+ = e + ^. This arrangement has the advantage that we can now 
write the ip-field in terms of four component spinors as : 
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The equations of motion become in matrix form 



(11) 
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-e_ Md -1 ^ 
-Md e_ -1 

1 -e+ -Md 

1 Md e+ J 



C2 



Mp 



C2 



(12) 



All the parameters are normalised by M. Inside the Q-ball the time dependent solution is : 






(13) 



where the up superscript stands for the first two components of the eigenvectors, while the 
down one indicates we take the two last components of the eigenvectors of the inner motion 
matrix. The time dependent solution outside the Q-ball is once more given by, 

4 

~Pj^ ' V J) jy- y-pj 



m = J2{Aj,Bj)e'^'-'^'^'u^Pe'P^' + {A*, B*)e-'^'+'^'^Hu'!°'^''e'P^'y 



(14) 



i=i 



where we this time use the eigenvectors of the outside motion matrix (Md = 0). The A's 
5's and C's are expension coefficients. The matching rules for these coefficients will give 
the reflection and transmission coefficients. 



III. INTERACTIONS OF Q-BALLS AND MATTER. 



We know from previous work 



a 



that we can use two equivalent methods to solve the 
problem of particle creation from Q-balls. These two methods allowed us to obtain the 
exact quantysed solutions describing fermions being produced by a Q-ball. The method 
we shall use now will consist in using the total solutions outside the evaporation range 
e ^ [—^ + Md; ^ + Md]. The solutions we constructed are superpositions of partial waves 
containing positive and negative movers. 

A. Interaction of Q ball and matter, massless case. 

We now study the diffusion processes on a Q ball. In the first part of the construction shall 
not consider any evaporation for Q balls. This simplification will allow us in the first place 
to build and study the diffusion of a fermion on a Q ball that we expect to be very similar 
as diffusion on a potential barrier. The simplest aproach will be to consider all possible 
transmition trough the Q ball, that is we might have a transmitted particle or anti-particle 
due to the energy shift (the ^ factor in the time dependent solution). For simplicity we 

p 



FIG. 1: General transmission through a Q ball. 

shall consider the construction where there is no incoming anti-particle. 

We start by using the solution outside the Q-Ball for the simple massles case we have: 

_(^Pi\ _( Cf e-^(^+^)*e*(^+^)^ 
\ip2 ) \ (Cl') V(^-^)*e^(^-1^)^ 

we have also a relation for coefficients on the right and on the left given by 



(15) 




(JR \ ^ / ^-2ikJ _ ^2^2ik,l _(y^Q-2ikJ ^ ^^Q2ikJ 

For the calculation of evaporation rate we considered that both components where the same 
particles what we do now is consider both componets to be different particles (doing so 



we can construct the normal field containing both particles and anti-particles), we shall 
consider as a first study : 

UJn UJn 

e+ — >0^e> 

2 - - 2 

2 - - 2 ' 

it means we have the same type of particles on the letf hand side moving both of them 
towards the Q-Ball. Since we we have chosen the parameters to be the incident amplitudes 
we can choose that there is no incident anti-particle, C2 = 0. Using the matrix relations we 
find, 

' ~ l-al ~ """ 

The transformation amplitude we are looking for is the C2 coefficient. As long as k^ is 
complex (see chap three), we can easily see that this amplitude decreases as the Q-Ball's 
size increases. Here we considered the range where e > ^ and the solution is, 

foo I Afe)e-*(^+^)*e*(^+1')* 
^^ 1^ 

foo /A(e)T,pe-(^+^)*e^(^+^)* 
^R= I dt\ 

J'f \ A*(e)T;^e*(^-1')*e^(^-^)* 

where the L and R subscripts indicate on which side of the Q-Ball we are. At this point 
the A expansion factor is considered to be continuous so it can be upgraded to an operator 
later on. It is not the only range where transformation occurs we also have the other range 
g < _ii|i^ with this time incident particles coming from the right hand side of Q-Ball, we 
have 

^RJ~l-a^y ^^^-2ikJ _ ^^g2^fcj ^2^kA _ c,2e-2*fe' j ^ ^ 

we have this time Cj^ = 0, leading to the same solution for the transmisson amplitudes, 

f = - = T 

J-pp 1 — Q;2 PP' 



even if the amplitudes are the same, wich is a effect of symmetry, we shall use tilde symbols 
to distinguish both contributions. In fact we would like to keep in mind that this second 
contribution comes from waves having negative energy. These two contributions will mix 
together to build the total solution. We have this time for solution, 

J-oo I 5'^(e)Tle^(^-^)V(^-^)^ 




^^ = y_ d^\ ,^ , .„. .. .0. I- (16) 

Let us take a look at the pieces on the left, 

r- — 
J —oo 

where mi = I and U2 = \ I . We change the varible e ^ — e in the integrals on the 

negative range to obtain. 



/*oo p 00 

^^ = / rfeMie-^(^+^)V(^+^)M(e) + / deuie'^'~^^'e-'^'-^^%p{-e)B{-e) 

+ r deu2e-'^'+^^'e~'^'+^^%p{-e)B*{-e). 

The next thing to do is change the variable e + ^ = e' and e — ^ = e' to obtain, 

^L = r deu,e-'^e-'^Aie' - ^) + T deu,e-%-' %p{-^ - e)B{-'^ - e) 
Jluo Z Jo I I 



00 



^Ox j^i,, / , ^0^ 



+ / deM2e-*^V^^%p(-e' + -f)5'^(-e 



UJi) 



2 ' ' 2 



Due to the energy shift comming from the Q ball we have two different bounds on the 
integration e G [uo, 00 [ and e G [0, 00 [. These two ranges will be of great interest, since each 
one can be identified to a different particle. Before we construct the quantised solution we 
need to do the same transformations to the solution on right-hand side. 

^^ = r deuie-'^^'+^^'e'^'+^^TppAie) + f' deu2e'^'-^^'S'-^^'T;pA\e) 

+ r"^ deu2e^^^-^-^^'e'^^-^-^>B\e), 
J —00 



after the variable changes and the sign change we obtain, 
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oo 



-i^it-i^'zrr II ^On ./ / ^0^ 



+ rfe«ie-'V-'%,(e'-7)A(e'-^). 

We see here that the B{—e' + ^) and A{e' — ^) are in the same integration range and 
therefor will corespond to the same particles. We now have both solutions one on each side 
of the Q-Ball, the last step of construction will be quantization. Quantization will be done 
identifying the parts coresponding to particles and the parts coresponding to anti-particles. 
It is easy to identify the particles and the anti-particles since the integration domain is well 
known and definied we have, 



for the solution on the left and 

for the solution on the right hand side of the Q-Ball. During the construction we used, 

B{k) = a{k)e{k) + b\k)e{-k). 

We keep for the moment the / and r subscript on the operators to identify the side of the 
Q ball they correspond. During this construction we did not normalise any solution. The 
normalisation of solutions is a simple task y, and will not be needed in the rest of the 
discussion. We just need to keep in mind that all the waves packets are normalised in the 
way that the sum of the reflected amplitude and the transmitted one is equal to one. 

Now that we have the total solution in terms of creation and anihilisation operators we 
can compute anything we want, by applying the proper operator valued observable. We 
can see that the energy of the transmited wave is the same as the energy of the incoming 
wave, while the energy of the reflected anti-particle is shifted by a factor of uo- This shift in 
energy, one component shifted downwards and the other one shifted upwards is the reason 
why we found finite evaporation ranges. This property is due to the Yukawa coupling inside 



the Q-Ball. We can choose a state where there is no incoming particles on the right so our 
transmission coeficients can be identified to refiection ones. In fact this construction can be 
used in the full energy range. The difficulty will be to seperate the phenomenons linked to 
diffusion and those linked to evaporation. 

B. Diffusion on a Q ball 

Now that we have the quantised solutions we can try to study diffusion on a Q ball. We 
need to use the total solution having no incident wave on the right-hand side of the Q ball. 
We write, 

+ l'^ rfee-*{e-%(^ + t)uM{^ + y ) 
+ r dee-\e-T,,{^ - e)u,al{^ - e)} 

the construction of the Bogolubov transformation linking the far past, the incident wave, to 
the far future reads : 

Cout = Tpp^Y - e)ai(e ^ y) + ^pp(y + Mi^ + y) 

We here distinguish two energy ranges, if e > ^ the aj(^ — e) operator becomes ar{e — ^), 
while if e < ^ it is the a/(e — ^) becomming aj(^ — e). As we said before we had to 
separate both regimes in order to see what type of diffusion we have. 

First we shall consider an infinite Q ball, so the right-hand side does not exist, and we 
have 

Cout = Tpp{^ - e)ai{e - y ) + ^pp(y + e)6j(e + y), 

(17) 

when we constructed the solution the incident particle had energy bigger than ujq. So in this 
range e — ^ is allways bigger than zero, and the a operator remains an anihilation operator. 
We have now in the final state (in the far future), 

< 0\c^c\0 >out= \Tpp\^ <p\p> . (18) 

10 



out 



This simply is a particle being reflected into an anti-particle to conserve helicity. To study 
the rest of the possibilities we shall use a small sized Qball. 

The first possibility stands for e > oj^jl we need to change the type of the second a 
operator, aj(^ — e) becomes a^[e — ^) : 



.^0 



^0^ 



.^0 



a;o^ 



Cout = Tpj,{— - t)ai[e - y) + Tpp[— + t)b\ (e + y 



,^0 



^0^ 



,^Q 



C^0^ 



Applying this to the vaccum state in the far future we have, 



out 



< 0|C^C|0 >o„i= \Tpp\^i < p\p >i + \Tpp\'^r < P\P >r- 



(19) 



FIG. 2: Diffusion on a Q ball for an incident particle having energy bigger than luq. 

The second possibility will be for a particle having energy lower than ujq/2 we have to 
change the type of the first a, a;(e — ^) becomes a[(^ — e). We write : 



,UJo 



UJQ. 



Cout = Tpp{— - e)al (y - e) + Tpp{— + e)bl (e + y J 



.^0 



,t^o 



.^0 



t^0^ 



+ Tppq - e)aiq -e)+ Tpp^ + e)blie + ^). 

This time we have four particles in the final state, it is the combination of evaporation and 
diffusion. 
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FIG. 3: Diffusion on a Q ball for an incident particle having energy smaller than uJo/2. 
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The third possibhty would be for incident fermions having their energy in the range 
e G [i^o/2\ cjq], we shall not study this range since we only considered incident fermions with 
enrgy bigger than uoq or smaller than uJq/2. But this range does not introduce any new or 
difficult construction. These construction can also be used to compute evaporation rate we 
just need to separate both constructions and then identify wich contribution is diffusion and 
wich one is evaporation. 

C. Diffusion on a Q ball massive case. 

Starting with the solution for a Q ball interacting with a fermion we have, 



iMo- 2 



CO 






Mo-^ 



as before this solution has two different integration ranges one for each type of particles. 
We shall now only consider the range where e > M/) + ^ in this range we only have diffusion 
and to not need to study the interaction of both, diffusion and evaporation. We have , 

the only reason why we used a very big Q ball is that we have less terms to deal with. 
With this solution we can as in the previous chapter construct a bogolubov transformation 
considering that in the far future only out going waves survies. We have this time : 

^ll^k + ^12^Pi + ^2l42 + ^22«P2 = Ccmt- (20) 

If we want Cout to have the same anti-commutation relations as a and h we need to normalise 
the operators with : 

(21) 

(22) 



^22 + ^21 



|n2r + kill 
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There is no obvious reasons why |ri2p = |^2iP and jrup = |r22p- These equahties will 
depend on the type of Q ball we study, they depend on the way the Q ball mixes the 
particles and anti-particle energies. We can use this c operator to compute the number of 
particles in the final state when the original one was the vacuum state, 

out < 0|c^c|0 >out = \r21\L < 0\aa^O >i„ +|rn|2„ < 0\bb^O >i„, 

= k2iP < p\p > +|^iiP < ap\ap > . (23) 

We have this time, when fermions are massive, a particle and an anti-particle in the final 
state. Those two refiected particules will be part of the state combining evaporation and 
diffusion. This combination of evaporation and diffusion is more complicated to study since 
we need to use the full solution and be carefull with the energy ranges and the type of 
operator linked to them. 

IV. CONCLUSIONS. 

As expected the interaction of Q balls and fermions can be separated into two different 
cases. The first one stands for interaction between the Q ball and fermions having their 
energy lying outside the evaporation rate. In this case we demonstrated that the interaction 
of Q balls and matter is nothing but standard diffusion. If we have an incident fermion 
we have a reflected and transmitted anti-fermion. The transmitted particles has its energy 
shifted by a uJo/2 factor due to interaction with the Q ball. 

The second case happens when the fermion interacting with the Q ball has its energy lying 
inside the evaporation range. In this case we have a superposition of two phenomenoms the 
flrst one is diffusion while the second one is evaporation. It seems that both phenomenoms 
to not interfear together and be considered separately. This fact is important in the way 
that it provides a new approach to compute Q ball evaporation rates. This new approach 
would be to study the diffusion of a particle on the Q ball's surface and flnd the two ranges 
where we have diffusion and both diffusion and evaporation. Then we just need to substract 
both amplitudes to obtain evaporation range. 

The interaction with massive fermions does not introduce any new facts. The calculations 
becomes more complicated but the separation into two ranges stands the same way, exept 
this time the evaporation range is different since it depends on the fermion mass. In this 

13 



case the other difference comes from the fact that we do not use hehcity conservation to find 
the refiected waves. 
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